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ABSTRACT 

We have analyzed galaxy and group-sized dark matter halos formed in a high 
resolution ACDM numerical N-body simulation in order to study the rotation of 
the triaxial figure, a property in principle independent of the angular momen- 
tum of the particles themselves. Such figure rotation may have observational 
consequences, such as triggering spiral structure in extended gas disks. The ori- 
entation of the major axis is compared at 5 late snapshots of the simulation. 
Halos with significant substructure or that appear otherwise disturbed are ex- 
cluded from the sample. We detect smooth figure rotation in 278 of the 317 halos 
in the sample. The pattern speeds follow a log normal distribution centred at 
Qp = 0.148 h km s~^ kpc~^ with a width of 0.83. These speeds are an order of 
magnitude smaller than required to explain the spiral structure of galaxies such 
as NGC 2915. The axis about which the figure rotates aligns very well with the 
halo minor axis, and also reasonably well with its angular momentum vector. 
The pattern speed is correlated with the halo spin parameter A, but shows no 
correlation with the halo mass. The halos with the highest pattern speeds show 
particularly strong alignment between their angular momentum vectors and their 
figure rotation axes. The figure rotation is coherent outside 0.12 r^ij.. The mea- 
sured pattern speed and degree of internal alignment of the figure rotation axis 
drops in the innermost region of the halo, which may be an artifact of the nu- 
merical force softening. The axis ratios show a weak tendency to become more 
spherical with time. 

Subject headings: dark matter — galaxies: evolution — galaxies: formation — 
galaxies: individual (NGC 2915) — galaxies: structure — galaxies: kinematics 
and dynamics 
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1. Introduction 

Although there have been many theoretical studies of the shapes of cosmological dark 
matter halos (e.g. Dubinski & Carlberg 1991; Warren et al. 1992; Cole & Lacey 1996; Jing 
& Suto 2002), there has been relatively little work done on how those figure shapes evolve 
with time, in particular, whether the orientation of a triaxial halo stays fixed, or whether the 
figure rotates. While the orientation of the halo can clearly change during a major merger, 
it is not known whether the orientation changes in between cataclysmic events. Absent any 
theoretical prediction one way or the other, it is usually assumed that the figure orientation 
of triaxial halos remain fixed when in isolation (e.g. Subramanian 1988; Johnston et al. 1999; 
Lee & Suto 2003) 

Early work at detecting figure rotation in simulated halos was done by Dubinski (1992) 
(hereafter D92). While examining the effect of tidal shear on halo shapes, he found that in 
all 14 of his 1-2 x 10^^ Mq halos, the direction of the major axis rotated uniformly around 
the minor axis. The period of rotation varied from halo to halo, and ranged from 4 Gyr at 
the fast end to 50 Gyr at the slow end, or equivalently the angular velocities ranged between 
0.1 and 1.6 km s^^ kpc~^. ^ It is difficult to draw statistics from this small sample size, 
especially since the initial conditions of this simulation were not drawn from cosmological 
models, but were performed in a small isolated box with the linear tidal field of the external 
matter prescriptively superimposed (Dubinski & Carlberg 1991). Given that the main result 
of D92 is that the tidal shear may have a significant impact on the shapes of halos, it is 
clearly important to do such studies using self-consistent cosmological initial conditions. 

Recent studies of figure rotation come from Bureau et al. (1999) (BFPM99) and Pfitzner 
(1999) (P99). P99 compared the orientation of cluster mass halos in two snapshots spaced 
500 Myr apart in an SCDM simulation (Q = 1, A = 0, /i = 0.5). He detected rotation 
of the major axis in ~ 5% of them, and argued that the true fraction with figure rotation 
is probably higher. BFPM99 presented one of these halos which was extracted from its 
cosmological surroundings and left to evolve in isolation for 5 Gyr. During that time, the 
major axis rotated around the minor axis uniformly at all radii at a rate of 60° Gyr~^, or 
about 1 km kpc~^. 

There may be observational consequences to a dark matter halo whose figure rotates. 
BFPM99 suggested that triaxial figure rotation is responsible for the spiral structure of the 
blue compact dwarf galaxy NGC 2915. Outside of the optical radius, NGC 2915 has a large 



^It may be more intuitive to think of angular velocity in units of radians Gyr ^ rather than the common 
unit of pattern speed, km kpc~^. Fortunately, the two units give almost identical numerical values. 
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H I disk extending to over 22 optical disk scale lengths (Meurer et al. 1996). The gas disk 

shows clear evidence of a bar, and a spiral pattern extending over the entire radial extent of 
the disk. BFPM99 argue that the observed gas surface density is too low for a bar or spiral 
structure to form by gravitational instability, and that there is no evidence of an interacting 
companion to trigger the pattern. They propose that the pattern may instead be triggered 
by a rotating triaxial halo. 

Bekki & Freeman (2002) followed this up with Smoothed Particle Hydrodynamics (SPH) 
simulations of a disk inside a halo whose figure rotates, and showed that a triaxial halo with a 
fiattening ofb/a = 0.8 and a pattern speed of 3.84 km s~^ kpc~^ could trigger spiral patterns 
in the disk. Masset & Bureau (2003) (hereafter MB03) found that in detail, the observations 
of NGC 2915 are better fit by increasing the disk mass by an order of magnitude (for example, 
if most of the hydrogen is molecular, e.g. Pfenniger, Combes, & Martinet 1994), but that a 
triaxial halo with b/a ^ 0.85 and a pattern speed of between 6.5 and 8.0 km kpc~^ also 
provides an acceptable fit. 

MB03 concluded that if the halo were undergoing sohd body rotation at this rate, its 
spin parameter would be A 0.157, which is extremely large (only 5 x 10~^ of all halos have 
spin parameters at least that large). However, this argument may be fiawed because the 
figure rotation is a pattern speed, not the speed of the individual particles which constitute 
the halo, and so it is in principle independent of the angular momentum; in some cases 
the figure may even rotate retrograde to the particle orbits (Freeman 1966). Schwarzschild 
(1982) discusses in detail the orbits inside elliptical galaxies with figure rotation. He finds 
that models can be constructed from box and X-tube orbits, which have no net streaming of 
particles with respect to the figure (though they have prograde streaming at small radius and 
retrograde streaming at large radius) , and so result in figures and particles with the same net 
rotation. He also constructs models that include prograde-streaming Z-tube orbits, which 
result in a figure that rotates slower than the particles. Stable retrograde Z-tube orbits also 
exist, but Schwarzschild (1982) did not attempt to include them in his models, so it may also 
be possible for the figure to rotate faster than the particles. While these results demonstrate 
the independence of the figure and particle rotation, it is not clear if they can be translated 
directly to dark matter halos. Dark matter halos may have different formation mechanisms 
and may be subject to different tidal forces than elliptical galaxies, and the different density 
profile may also have a large effect on the viable orbital families (Gerhard & Binney 1985). 

If there are observational consequences to dark halo figure rotation, such as those found 
by Bekki & Freeman (2002), they can be used as a direct method to distinguish between dark 
matter and models such as MOdified Newtonian Dynamics (MOND) that propose to change 
the strength of the force of gravity (Milgrom 1983; Sanders & McGaugh 2002). Many of the 
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traditional methods of deducing dark matter cannot easily distinguish between the presence 
of a roughly spherical dark matter halo and a modified force or inertia law. However, a major 
difference between dark matter and MOND is that dark matter is dynamical, and so tests 
that detect the presence of dark matter in motion are an effective tool to discriminate between 
these possibilities. Among the tests that can make this distinction are the ellipticities of dark 
matter halos as measured using X-ray isophotes, the Sunyaev-Zeldovich effect, and weak 
lensing (Buote et al. 2002; Lee & Suto 2003, 2004; Hoekstra et al. 2004); the presence of 
bars with parameters consistent with being stimulated by their angular momentum exchange 
with the halo (Athanassoula 2002; Valenzuela & Klypin 2003); and spatial offset between 
the baryons and the mass in infalling substructure measured using weak lensing (Clowe, 
Gonzalez, & Markevitch 2004). Rotation of the halo figure requires that dark matter is 
dynamic, and therefore observable structure triggered by figure rotation potentially provides 
another test of the dark matter paradigm. 

In this paper, we use cosmological simulations to determine how the figures of ACDM 
halos rotate. The organization of the paper is as follows. Section 2 presents the cosmological 
simulations. Section 3 describes the method used to calculate the figure rotations, which are 
presented in Section 4. Finally we discuss our conclusions in Section 5. 



2. The Simulations 

The halos are drawn from a large high resolution cosmological N-body simulation per- 
formed using the GADGET2 code (Springel et al. 2001). We adopt a "concordance" cos- 
mology (e.g. Spergel et al. 2003) with = 0.3, = 0.7, (^bar = 0.045, h = 0.7, and 
(Tg = 0.9. The only effect of r2bar is on the initial power spectrum, since no baryonic physics 
is included in the simulation. The simulation contains 512^ = 134, 217, 728 particles in a 
periodic volume 50 Mpc comoving on a side. All results are scaled into /i-independent 
units when possible. The full hst of parameters is given in Table 1. 

A friends-of-fricnds algorithm is used to identify halos (Davis et al. 1985). We use the 
standard linking length of 

b = 0.2f^-^/^ (1) 

where n = N/V is the global number density. 

Measuring the figure rotation requires comparing the same halo at different times during 
the simulation. We analyze snapshots of the simulation at lookback times of approximately 
1000, 500, 300, and 100 h^^ Myr with respect to the z = {) snapshot. The scale factor a of 
each snapshot, along with its corresponding redshift and lookback time, is listed in Table 2. 
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Table 1. Parameters of the cosmological simulation. 



Parameter Value 



N 5123 

Box size {h~^ Mpc comoving) 50 

Particle mass (10^ h-^ Mq) 7.757 
Force softening length {h~^ kpc) 5 
Hubble parameter h {Ho = 100 h km s"^ Mpc"^) 0.7 

Qo 0.3 

Qa 0.7 

(78 0.9 

Jlbar 0.045 



Table 2. Snapshots used to calculate figure rotations. 



Snapshot Name 


Scale Factor 


Redshift 


Lookback Time 






(--) 


Mat) 


b090 


0.89 


0.1236 


1108 


b096 


0.95 


0.0526 


496 


b098 


0.97 


0.0309 


296 


blOO 


0.99 


0.0101 


98 


bl02 


1.0 


0.0 
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3. Methodology 
3.1. Introduction 

The basic method is to identify individual halos in the final z — snapshot of the 
simulation, to find their respective progenitors in slightly earlier snapshots, and to measure 
the rotation of the major axes through their common plane as a function of time. 

Precisely determining the direction of the axes is crucial and difficult. When merely cal- 
culating axial ratios or internal alignment, errors on the order of a few degrees are tolerable. 
However, if a pattern speed of 1 km kpc~^, as observed in the halo of BFPM99, is typical, 
then a typical halo will only rotate by 4° in between the penultimate and final snapshots of 
the simulation. Therefore, the axes of a halo must be determined more precisely than this 
in order for the figure rotation to be detectable. In fact, we should strive for even smaller 
errors to see if slower-rotating halos exist. It would have been difficult for P99 to detect halos 
rotating much slower than the halo presented in BFPM99; although the error varies from 
halo to halo (for reasons discussed in section 3.3), Figure 5.23 of P99 shows that most of his 
halos had orientation errors of between 8° and 15°, corresponding to a minimum resolvable 
ffgure rotation of ~ 0.6 km s~^ kpc~^ for a 2a detection in snapshots spaced 500 Myr apart. 

A major difficulty in determining the principal axes so precisely is substructure. The 
orientation of a mass distribution is usually found by calculating the moment of inertia 
tensor lij = J2k''^ki^k,ii^k,j, and then diagonalizing I^j to find the principal axes. However, 
this procedure weights particles by r^. Therefore, substructure near the edge of the halo (or 
of the subregion of the halo used to calculate the shape) can exert a large infiuence on the 
shape of nearly spherical halos, especially if a particular subhalo is part of the calculation 
in one snapshot but not in another, such as when it has just fallen in. This is particularly 
problematic because subhalos are preferentially found at large radii (Ghigna et al. 2000; 
De Lucia et al. 2004; Gill et al. 2004; Gao et al. 2004). Moving substructures can also 
induce a false measurement of figure rotation due to their motion within the main halo at 
approximately the circular velocity. 

To mitigate this, we firstly use particles in a spherical region of radius 0.6 r^ii. surrounding 
the center of the halo, rather than picking the particles from a density dependent ellipsoid as 
in Warren et al. (1992) or Jing & Suto (2002). We find that those methods allow substructure 
at one particular radius to infiuence the overall shape of the ellipsoid from which particles 
are chosen for the remainder of the calculation, and therefore bias the results even when 
other measures are adopted to minimize their effect. The choice of a spherical region biases 
the derived axis ratios toward spherical values, but does not affect the orientation. Secondly, 
the particles are weighted by 1/r^ so that each mass unit contributes equally regardless of 
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radius (Gerhard 1983). Both D92 and P99 take similar approaches, but using radii based 
on elhpsoidal shells. Therefore, we base our analysis on the principal axes of the reduced 
inertia tensor 



In the majority of halos, the substructure is a small fraction of the total mass, usually less 
than 5% of the total mass within 60% of the virial radius (De Lucia et al. 2004, Figure 8), 
so its effect is much reduced. There are still some halos which have not yet relaxed from a 
recent major merger, in which case the "substructure" constitutes a significant fraction of 
the halo. To find these cases, the reduced inertia tensor is separately calculated enclosing 
spheres of radius 0.6, 0.4, 0.25, 0.12, and 0.06 times the virial radius to look for deviations 
as a function of radius (see section 3.4 for details). These radii are always with respect to 
the ^ = value of Tyir- 

We find that only halos with at least 4 x lO"' particles, or masses of at least ~ 3 x 
10^^ h'^ Mq have sufficient resolution for the orientation of the principal axis to be deter- 
mined at sufficient precision (see Section 3.3). There are 1432 halos in the z = snapshot 
satisfying this criterion, with masses extending up to 2.8 x 10^^ Mq. 



To match up the halos at 2; = with their earlier counterparts, we use the individual 
particle numbers provided by GADGET which are invariant from snapshot to snapshot, and 
find which halo each particle belongs to in each snapshot. The progenitor of each z = 
halo in a given z > snapshot is the halo that contributes > 90% of the final halo mass. 
Sometimes no such halo exists; in these cases, the halo has only just formed or underwent 
a major merger and so is not useful for our purposes. Figure 1 shows a histogram of the 
fraction of the final halo mass that comes from the b096 {z 0.05) halo which contributes 
the most mass. There are also some cases where two nearby objects are identified as a single 
halo in an earlier snapshot but as distinct objects in the final snapshot. We therefore impose 
the additional constraint that the mass contributed to the final halo must also be > 90% of 
the progenitor's mass. In the longer time between the earliest snapshot b090 and the final 
snapshot bl02, a halo typically accretes a greater fraction of its mass, and so a more liberal 
cut of 85% is used for this snapshot (see the dashed histogram in Figure 1). 492 of the halos 
that satisfied the mass cut did not have a progenitor which satisfied these criteria in at least 
one of the 2; > snapshots and so were eliminated from the analysis, leaving a sample of 940 
matched halos. 




(2) 



3.2. Halo matching 
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Fig. 1. — Histogram of the fraction of the final mass that comes from the b096 {z ~ 0.05, 
sohd fine) and b090 {z ~ 0.12, dashed fine) halo which contributes the most mass. 
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3.3. Error in ctxis orientation 

There are two sources of errors that enter into the determination of the axes: how well 
the principal axes of the particle distribution can be determined, and whether that particle 
distribution has a smooth triaxial figure. Here we estimate the error assuming that it is 
not biased by substructure. The hales for which this assumption does not hold will become 
apparent later in the calculation. 

For a smooth triaxial ellipsoid represented by particles, the error is a function of N 
and of the intrinsic shape: as the axis ratio b/a or c/b approaches unity, the axes become 
degenerate. To quantify this, we have performed a bootstrap analysis of the particles in a 
sphere of radius 0.6 Tyir of each z — halo (Heyl et al. 1994). If the sphere contains N 
particles then we resample the structure by randomly selecting N particles from that set 
allowing for duplication and determine the axes from this bootstrap set. We do this 100 
times for each halo. The dispersion of these estimates around the calculated axis is taken 
formally as the "Icr" angular error, and is labelled 6'boot- 

As expected, the two important parameters are N and the axis ratio. We focus here 
on the major axis, for which the important axis ratio is b/a. The top panels of Figures 2 
and 3 show the dependence of the bootstrap error on and b/a respectively for all halos 
with M > 10^^ h^^ Mq. The solid lines are empirical fits: 



^err,Af - — (3) 



and 



0err,b/a = O-OO^J^^- (4) 

The form of equation (3) is not surprising; if a smooth halo was randomly sampled, we would 
expect the errors to be Poissonian with an N^^^"^ dependence. However, the cosmological 
halos are not randomly sampled. Individual particles "know" where the other particles are, 
because they have acquired their positions by reacting in the potential defined by those other 
particles. Therefore, the errors may be less than expected from a randomly sampled halo. 
To test this, we construct a series of smooth prolate NFW halos (Navarro, Frenk, & White 
1996) with b/a axis ratios ranging from 0.5 to 0.9, randomly sampled with between 3 x 10^ 
and 3 x 10^ particles, and perform the bootstrap analysis identically for each of these halos 
as for the cosmological halos. Because the method for calculating axis ratios outlined in 
Section 3.1 biases axis ratios toward spherical, the recovered b/a of these randomly sampled 
halos is larger than the input value, and ranges from 0.65 to 0.95. The errors for these 
randomly sampled smooth halos are shown as asterisks in Figures 2 and 3. 
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Fig. 2. — Angular bootstrap error ^boot a function of the number of particles within the 
central 0.6 r^i,. of each halo. Points are the cosmological halos, and asterisks are randomly 
sampled smooth NFW halos. (Top): Angular error 6'boot- The solid line is the fit 6'err,Ar 
from equation (3). (Middle): Ratio between the angular error and the error expected for 
the halo given its axis ratio b/a, i.e. 6'boot /6'err,Va- The sohd line is 6'err,Ar from equation (3) 
renormalized to the typical error of 0.02 radians. (Bottom): Ratio between the angular error 
and the analytic estimate ^err from equation (5). 
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Fig. 3. — Angular bootstrap error ^boot a function of the axis ratio h/a of each halo. 
Points are the cosmological halos, and asterisks are randomly sampled smooth NFW halos. 
(Top): Angular error 6'boot- The solid line is the fit Oerr,b/a from equation (4). (Middle): 
Ratio between the angular error and the error expected for the halo given the number of 
particles N, i.e. O\^oot/0cTT,N- The solid line is 9eTT,b/a from equation (4) renormalized to the 
typical error of 0.02 radians. (Bottom): Ratio between the angular error and the analytic 
estimate 6'err from equation (5). 
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The top panel of Figure 2 shows a rise in the dispersion of the error for N < 4000, with 
many halos having errors greater than the 0.1 radians necessary to detect the figure rotation 
of the halo presented in BFPM99. Therefore, we only use halos with N > 4000. 

The bootstrap error appears to be completely determined by N and b/a. The residuals 
of ^boot with respect to ^err,Ar are due to 9err,b/a and vice versa. This is shown in the middle 
panels of Figures 2 and 3. In the middle panel of Figure 2 we have divided out the dependence 
of ^boot on the axis ratio, making apparent an extremely tight relation between the residual 
and N, while in the middle panel of Figure 3 we have divided out the dependence of ^^boot 
on A^, showing the equally tight relation between the residual and b/a. It is apparent from 
comparing the points and asterisks that the errors in the cosmological halos are slightly 
smaller than for randomly sampled smooth halos. 

Combining equations (3) and (4), and noting that the typical error is ^boot ~ 0.02 radi- 
ans, we find the bootstrap error is well fit by 



The bottom panels of Figures 2 and 3 show the residual ratio between the bootstrap error 
^boot and the analytic estimate ^^crr- The vast majority of points lie between 0.8 and 1.0, 
indicating that ^^crr overestimates the error by ~ 10%. Equation (5) breaks down as b/a 
approaches unity; these halos are nearly oblate and so do not have well-defined major axes. 
It also becomes inaccurate at very low b/a due to low-mass poorly-resolved halos. Even in 
these cases, the error estimate is conservative, but to be safe we have eliminated axes with 
b/a < 0.35 or b/a > 0.95 from the subsequent analysis, regardless of the nominal error. The 
randomly-sampled smooth halos follow equation (5) extremely well, so the non-Poissonianity 
of the sampling in simulated halos reduces the errors by 10%. 

Calculating the bootstraps is computationally expensive, so equation (5) is used for the 
error in all further computation. Because this estimate is expected to be correct for smooth 
ellipsoids, cases where the error is anomalous are indications of residual substructure. 



Ideally one would fit the figure rotation by comparing the orientation of the axis at each 
snapshot to that of a unit vector rotating uniformly along a great circle, and minimize the 
to find the best fit uniform great circle trajectory. However, this requires minimizing a 
non-linear function in a four-dimensional parameter space, a non-trivial task. We adopt the 
simpler and numerically more robust method of solving for the plane z — ax + by that fits the 




(5) 



2VNl-b/a' 



3.4. Figure rotation 
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Fig. 4. — Diagram that demonstrates how we fit a plane to the major axis measurements at 
all snapshots (thick lines) and then find the increase of phase as a function of time. The 
figure rotation axis is perpendicular to the best fit plane, and defined such that (p increases 
around it counter-clockwise with time. 
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major axis measurements of the halo best at all timesteps, assuming the error is negligible. 
The change of the phase of the axes in this plane as a function of time are then fit by linear 
regression. A schematic diagram of this process is shown in Figure 4. 

The degree to which the axes are consistent with lying in a plane is checked by calculating 
the out-of- plane y^: 



where v is the number of degrees of freedom and A^j is the minimum angular distance 
between the major axis at timestep i and the great circle defined by the best fit plane. 

Because the axes have reflection symmetry, it is impossible to measure a change in 
phase of more than 7r/2. The phases are adjusted by units of tt such that the difference in 
phase between adjacent snapshots is always less than 7r/2. If the figure were truly rotating 
by 90° or more in between the snapshots, it would be impossible to accurately measure 
this rotation since the angular frequency would be larger than the Nyquist frequency of our 
sampling rate. Any faster pattern speeds would be aliased to lower angular velocities, with 
an aliased angular velocity of f^Nyq — iS^p — ^Nyq), where VL^ is the intrinsic angular velocity 
of the pattern and f^Nyq is the Nyquist frequency. For snapshots spaced 500 h"^ Myr apart, 
the maximum time between the snapshots we analyze, the maximum detectable angular 
velocity is 3.8 h km s~^ kpc~^. We do not expect the figure to change so dramatically as we 
have excluded major mergers. However, this can be checked 'post facto by checking whether 
the distribution of measured angular velocities extends up to the Nyquist frequency; if so, 
then there are likely even more rapidly rotating figures whose angular frequency is aliased 
into the detectable range, fooling us into thinking they are rotating slower. If the measured 
distribution does not extend to the Nyquist frequency, then it is unlikely that there are any 
figures rotating too rapidly to be detected (see Section 4) . 

The best fit linear relation for the phase as a function of time is found by linear re- 
gression. Because the component of an isotropic angular error projected onto a plane is 
half of the isotropic error, we divide the error of equation (5) by two before we perform the 
regression. The slope of the hnear fit gives the pattern speed Qp of the figure rotation. The 
error is the la limit on the slope. 

Once we have calculated the pattern speed for each halo, we can eliminate the cases 
where substructure has severely impacted the results. In these cases, the signal is dominated 
by a large subhalo at a particular radius, so the derived pattern speed will be significantly 
different when the sphere is large enough to include the subhalo from when the subhalo is 
outside the sphere. We have calculated the pattern speed using enclosing spheres of radius 
0.6, 0.4, 0.25, 0.12, and 0.06 of the virial radius. The fraction of mass in subhalos can be 




(6) 
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Fig. 5. — Log-weighted projected density of 4 halos with a range of subhalo fractions fs- 
The subhalo fractions are 0.166 (top-left), 0.065 (top-right), 0.045 (bottom- left), and 0.016 
(bottom-right). Axes are in units of kpc from the halo center. All halos have masses in 
the range 2 - 3 x 10^^ h'^ M©. 
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estimated via the change in the pattern speed Qp at adjacent radii. Because the reduced 
inertia tensor is mass-weighted, the figure rotation of a sphere with a smooth component 
rotating at fip^smooth plus a subhalo containing a fraction fg of the total mass moving at the 
circular velocity Vc at radius R Tyir is approximately 

Qp !V {1 — /s)ilp,smooth + fs-^ , (7) 

where the difference due to the presence of the subhalo is 



lGM{< R 



If the density profile is roughly isothermal, the enclosed mass is 

M(< R rvir) = ^7rA,p,Hti? (9) 

Solving equations (8) and (9) gives an expression for the fraction of the mass in substructure 
given a jump of Aflp in the measured pattern speed when crossing radius R ryir' 

f. = ^ia^. (10) 

y |7rGAcPcrit 

The term in the square root is equal to 0.72 h km s~^ kpc~^. For each halo, we compute the 
value of fs due to the jump Aflp between each set of adjacent radii, i.e. for R — 0.4, 0.25, 0.12, 
and 0.06, and adopt the largest value of fg as the substructure fraction of the halo. Figure 5 
shows log-weighted projected densities of 4 halos in the mass range 2 - 3 x 10^^ Mq with 
a variety of values of fs, ranging from 0.166 at the top-left to 0.016 at the bottom right. 
After examining a number of halos spanning a range of fs, we adopt a cutoff of fs < 0.05 
for undisturbed halos. This eliminates 289 of the 940 halos, leaving 651 undisturbed halos. 

A further 158 halos were eliminated because the angular error approached 7i/2 in at 
least one of the snapshots. This includes the halos with b/a < 0.35 or 6/a > 0.95 discussed 
in Section 3.3. We also eliminate cases where the reduced from the linear fit of phase 
versus time indicates that the intrinsic error of the direction determination is much lower 
than suggested by equation (5), indicating that the model of the halo as a smooth ellipsoid 
is violated (10 halos with xl < 0.1), and those cases where the phase does not evolve linearly 
with time (134 halos with xl > 10). Finally, we ehminate halos where the axes do not he 
on a common plane, i.e. the 32 halos where xlop ^0. Therefore, the final sample consists 
of 317 halos. 



A sample halo is shown in the first five panels of Figure 6. It was chosen randomly 
from the halos with relatively low errors and typical pattern speeds. It has a mass of 
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Fig. 6. — (Upper-left five panels): Projection onto the best fit plane of the inner 0.6 Tyir 
of a sample halo at the five snapshots we analyze. Axes are in units of kpc from the 
halo center. From left to right, top to bottom, the snapshots are at 1108, 496, 296, 98, and 
Myr before z = 0. The solid line is the major axis, which rotates counterclockwise by 
20° from beginning to end. (Bottom-right): Phase of the major axis in the rotational plane 
of the sample halo. The zero point is arbitrary, but identical in all snapshots. The solid line 
is the hnear fit, with a slope of 0.33 h radians Gyr"^. 
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1.9 X 10^^ Mq, and a pattern that rotates at 0.32 ± 0.01 h km kpc~^. It has a spin 
parameter A = 0.047, and axis ratios of b/a = 0.86 and c/a = 0.77 ai z = 0. The derived 
substructure fraction is fg = 0.045, and the out-of-plane xlop = 8.5. The sohd hne shows 
the measured major axis in each snapshot, which rotates counterclockwise in this projection. 
The phase of its figure rotation as a function of time is shown in the bottom-right panel 
of Figure 6. The zero point is arbitrary, but is consistent from snapshot to snapshot. The 
hnear fit is also shown, which has a reduced of 2.9. 



Figure 7 shows the measured figure rotation speeds for all of the halos in the sample, 
as a function of their error. Halos with measured pattern speeds less than twice as large 
as the estimated error (the dashed line) are taken as non-detections. 278 of the 317 halos 
have detected figure rotation. A histogram of the pattern speeds is presented in Figure 8, 
expressed in logJlp. The thin histogram contains all halos with 2a detections, while the 
thick histogram contains those with the smallest errors, less than 0.01 h km s~^ kpc~^. The 
largest upper limit due to a non-detection in the main sample is flp < 0.126 (logJlp < —0.9), 
so the thin histogram is incomplete below this level, while the low error sample contains only 
one upper limit, Qp < 0.015 {logQp < —1.8), so the thick histogram is complete down to this 
level. The dashed curves are Gaussian fits to the histograms. The fit to the thin histogram, 
which has the largest sample size but is incomplete at low flp, peaks at log flp = —0.80 and 
has a standard deviation of 0.29, while the thick curve, which contains fewer halos but is less 
biased toward large values of flp, peaks at log flp — —0.84 and has a standard deviation of 
0.34. We give more weight to the thick histogram, whose points all have very small errors, 
and propose that the true distribution peaks at log flp — —0.83 with a standard deviation of 
0.36. Expressed as a log normal distribution, the probability is 



This fit is shown in Figure 9, compared to the fractional distribution of halos in the full 
(thin) and low error (thick) samples, and encompasses both the large number of halos with 
low flp seen when the errors are sufficiently small, and the tail at high flp seen when the 
sample size is sufficiently large. 

For comparison, the halo in BFPM99 has a pattern speed of 2 /i km s"^ kpc~^. This 
lies slightly above the top end of our distribution; the maximum pattern speed in our sample 
is 1.01 h km s"^ kpc~^. Based on the log normal fit of equation (11), we estimate the 



4. 



Results 




(11) 



where flp^ = iq-o-ss ^ o.l48 h km s'^ kpc"^ and the natural width a = 0.36 In 10 = 0.83. 
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Fig. 7. — Measured pattern speed of the rotation of the major axis. The x-axis is the error 
in the pattern speed. The sohd hne shows where the measure pattern speed is equal to the 
estimated error, while the dashed line is the 2a limit. The horizontal dotted line shows the 
Nyquist limit. 
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Fig. 8. — Histogram of the pattern speeds of the figure rotation, expressed in log Vtp. The 
thin histogram contains all halos that have 2a detections of figure rotation, i.e. those above 
the dashed fine of Figure 7, and is incomplete at < 0.126 h km kpc~^ or equiva- 
lently logQp < —0.9. The thick histogram contains only those halos with errors less than 
0.01 h km s~^ kpc~^, and is incomplete at < 0.015 h km s~^ kpc~^ or logQp < —1.8. 
The dashed curves are Gaussian fits to the histograms. The vertical dotted line shows the 
Nyquist limit. 



c 

o 



o 
o 




0.4 0.6 0. 

fip {h km s"'' kpc"^) 



1 .2 



Fig. 9. — Fractional histogram of pattern speeds of the figure rotation, expressed hnearly 
in Q,p. The thin histogram contains all halos that have 2a detections of figure rotation 
(incomplete at Jlp < 0.126 h km s~^ kpc~^), while the thick histogram contains halos with 
errors less than 0.01 h km s~^ kpc~^and is complete down to Q,p — 0.015 h km s~^ kpc~^. 
The dashed curve is the log normal fit given by equation (11). 
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fraction of halos with Qp > 2 h km kpc~^ to be ~ 10~^. Therefore, this halo is unusual, 
but it is not unreasonable to find a halo with such a pattern speed in a large simulation. 
Given the size of the errors in P99, and that he found very few halos with figure rotation, it 
should not be surprising that P99 could only detect pattern speeds at the upper end of the 
overall distribution. The different adopted cosmologies may also influence the results (note, 
however, that this comparison is performed in /i-independent units). Our results are also 
mostly consistent with D92, who finds pattern speeds of between 0.1 and 1.6 km s"^ kpc~^ 
in a sample of 14 halos. We have trouble reproducing the most rapidly rotating halos in 
D92, but this may be a product of the heuristic initial conditions in D92 compared to the 
cosmological initial conditions we use. 

In order to account for the spiral structure in NGC 2915, a triaxial figure would need to 
rotate at 7± 1 km s~^ kpc~^ (MB03). This is almost an order of magnitude faster than the 
fastest of the halos in our sample, and the log normal fit from equation (11) suggests that 
the fraction of halos with Qp >6 km s~^ kpc~^ is 5 x 10~^. Therefore, the figure rotation of 
undisturbed ACDM halos can not explain the spiral structure of NGC 2915. SPH simulations 
of gas disks inside triaxial halos with pattern speeds of 0.77 km s^^ kpc^^, comparable to the 
fastest pattern speeds in our sample, show very weak if any enhancement of spiral structure 
compared to a static halo (Bekki & Freeman 2002, Figure 2f). Therefore, it is unlikely that 
triaxial figure rotation can be detected in extended gas disks. 

The dotted line in both Figures 7 and 8 shows the Nyquist frequency of 3.8 h km s^^ kpc^ 
If the measured distribution of pattern speeds extended up to the Nyquist frequency, the 
intrinsic distribution would likely extend above the Nyquist frequency, and the results would 
be affected by frequency aliasing. However, the measured distribution does not approach 
the Nyquist frequency. Therefore, any halo whose figure rotation is aliased would need to 
be wildly anomalous, with a pattern speed many times faster than any other halo in our 
sample. We consider this unhkely. Figure 10 shows the pattern speeds as a function of the 
error, as in Figure 7, except that it only uses snapshots b096 through bl02, so the maxi- 
mum time between snapshots is 200 h^^ Myr, and the corresponding Nyquist frequency is 
7.6 h km s~^ kpc^^, shown again as the dotted line. The top of the distribution does not 
change between Figures 7 and 10, demonstrating that the results are not affected by aliasing. 

We investigate how the figure rotation axis relates to two other important axes. Both 
D92 and P99 claim that the major axis rotates around the minor axis. The direction cosine 
between the rotation axis and the minor axis is plotted both as a function of the pattern 
speed and as a histogram in Figure 11. We confirm that the rotation axis coincides very 
strongly with the minor axis. Due to the tendency for the halos to be prolate, the minor and 
intermediate axes tend to degenerate, accounting for some (but not all) of the misaligned 
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Fig. 10. — As in Figure 7, but only including snapshots b096 through bl02. 
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Fig. 11. — (Top): Histogram of the direction cosine between the minor axis of the halo 
and the figure rotation axis, for halos with 2a detections of figure rotation. Because the 
minor axis has reflection symmetry, this is always positive. (Bottom): Direction cosine as a 
function of the magnitude of figure rotation. 
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tail. 

The rotation axis is compared to the angular momentum vector of the halo in Figure 12. 
Because the angular momentum is usually well aligned with the minor axis of halos (Warren 
et al. 1992), it is no surprise that the rotation axis is also well aligned with the angular 
momentum vector. Because the alignment between the minor axis and angular momentum 
of halos is not perfect (Warren et al. (1992), Bailin & Steinmetz, in preparation), some of 
the halos with perfectly aligned figure rotation and minor axes have less perfect alignment 
between the figure rotation axis and the angular momentum vector, seen as the bump in 
Figure 12 that extends down to a direction cosine of 0.6; the tail of the distribution extends 
all the way to ant i- alignment. The alignment is also plotted as a function of the pattern 
speed in the lower panel of Figure 12. There is no trend for the halos with slow figure 
rotation, but all but two of the halos with Qp > 0.4 h km s~^ kpc~^ have figure rotation 
axes and angular momentum vectors that are well aligned, with a direction cosine of 0.65 or 
higher. 

We have attempted to see if the pattern speed is correlated with other halo properties, 
in particular its mass and its angular momentum. Figure 13 shows the pattern speed of the 
figure rotation versus the halo mass. Error bars are la errors, with 2a upper limits plotted 
for halos which lie below the dashed line of Figure 7. There is no apparent correlation 
between the halo mass and its pattern speed. 

Figure 14 shows the pattern speed versus the spin parameter A, where 



J\E 



1/2 



(Peebles 1969). We use the computationally simpler A' as an estimate for A, where 

A' = -J- (13) 

y/2MVR 

(Bullock et al. 2001). There is a tendency for halos with fast figure rotation to have large 
spin parameters; in particular, all but one of the halos with > 0.4 h km s~^ kpc"^ have 
A > 0.024. These are the same halos which are shown to have particularly well-aligned 
rotation axes and angular momentum vectors in Figure 12. We have calculated the median 
value of Vlp including the upper limits for bins of width AA = 0.01. The median rises steadily 
from 0.12 h km s^^ kpc'^ for A < 0.02 to 0.44 h km s^^ kpc"^ for A > 0.06. Note that P99 
only detected figure rotation in halos with A > 0.05 (see his Figure 5.24). 

The degree of alignment between the figure rotation axis and the angular momentum 
vector may depend on the angular momentum content of the halo. Figure 15 shows how 
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Fig. 12. — (Top): Histogram of the direction cosine between the angular momentum vector 
of the halo and the figure rotation axis, for halos with 2a detections of figure rotation. 
(Bottom): Direction cosine as a function of the magnitude of figure rotation. 
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Fig. 13. — Pattern speed of the figure rotation versus the mass of the halo. Error bars are 
1(7 errors for halos with at least 2a detections of figure rotation. The upper limits are the 
halos with a measured Qp < 2(7, and are plotted at the 2(7 limit. 
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Fig. 14. — Pattern speed of the figure rotation versus the spin parameter of the halo. Error 
bars and upper limits are as in Figure 13. 
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Fig. 15. — Direction cosine between the angular momentum vector of the halo and the figure 
rotation axis versus the spin parameter A. 
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this alignment depends on the spin parameter A. There is no trend for A < 0.05, but the 
halos with A > 0.05 show particularly good alignment. This is a natural consequence of 
the tendency for halos with rapid figure rotation to have well aligned figure rotation and 
angular momentum axes (Figure 12), and the correlation between figure pattern speed and 
spin parameter A (Figure 14). 

Figures 16 and 17 show how the figure rotation changes with radius. Figure 16 shows 
how the pattern speeds at different radii are related, while Figure 17 shows the alignment of 
the figure rotation axes between radii. Each panel includes only the halos that have at least 
4000 particles within the inner radius, pass all of the tests of Section 3.4 for both radii, and 
have 2(7 measurements of figure rotation at both radii. Due to the smaller number of particles 
in spheres of smaller radii, there are progressively fewer halos with good measurements at 
smaller radii. The top panels show that the figure rotation in the outer regions of the halo is 
very coherent. To some degree, this is by construction; the test for substructure is equivalent 
to a cut in AQp between adjacent radii. However, gradual drifts of flp and changes in the 
figure rotation axis with radius are still possible. The lower panels of Figures 16 and 17 
show that this indeed happens. In particular, while the figure rotation within 0.12 rvir and 
0.6 Tvir are strongly correlated, the pattern speeds within 0.12 r^ij. are slightly smaller and 
the alignment of the rotation axes is not quite as strong. The bottom panels show that in 
the innermost regions, within 0.06 ryir, the pattern speeds are significantly smaller than for 
the halo as a whole, particularly for those halos with high pattern speeds, and more than 
half of the halos show no alignment between the figure rotation axes. 

We examine three possible explanations for these trends with radius. First, it may be 
that the halos with high pattern speeds are still affected by residual substructure in the outer 
regions. However, the gradual dechne for all halos seen as the radius shrinks suggests that 
the mechanism responsible for the difference affects all halos equally and gradually, rather 
than affecting a few halos at a specific radius. Another piece of evidence that argues against 
this explanation is that the halos with the highest measured pattern speeds do not have 
preferentially high values of fg] they have values evenly spread between and the cutoff of 
0.05. A second possibility is that the figure rotation, although steady on timescales of 1 Gyr, 
may be fundamentally a transitory feature caused by a tidal encounter or the most recent 
major merger. The inner region of the halo has a shorter dynamical time, and therefore the 
effects of such a disturbance will be erased faster in the inner regions than the outer regions. 
This is consistent with the gradual decrease in pattern speed with radius and the decrease 
in ahgnment. However, the halo of BFPM99 has fast figure rotation (faster than any of our 
halos), and yet shows steady figure rotation at all radii for 5 Gyr. We propose instead that 
the effects of force softening are becoming important at the smaller radii. The radius of 
the innermost sphere for the halos plotted in the bottom panels of Figures 16 and 17 range 
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Fig. 16. — Pattern speed of figure rotation Qp at 0.4, 0.25, 0.12, and 0.06 of the virial radius 
Tvij. (top to bottom) as a function of the pattern speed at 0.6 Tyir- Only halos where both 
radii in the comparison contain at least 4000 particles, pass all of the tests of Section 3.4, 
and have 2a detections of figure rotation arc included. All units are h km s~^ kpc~^. The 
solid line corresponds to equal pattern speeds. 
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Fig. 17. — Histogram of the direction cosine of the ahgnment between the figure rotation 
axis of the halo as a whole (r < 0.6 r^i^) and of just the inner 0.4, 0.25, 0.12, and 0.06 of the 
virial radius Tvir (top to bottom). Only halos where both radii in the comparison contain 
at least 4000 particles, pass all of the tests of Section 3.4, and have 2a detections of figure 
rotation are included. 
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from 3-5 force softening lengths, where the effects of the gravitational softening can still be 
important (Power et al. 2003). The weaker gravitational force results in a more spherical 
potential, consistent with the weaker figure rotation and lack of alignment. 

We have calculated the rate of change of the h/a and c/a axis ratios over the 5 snapshots 
using linear regression. The evolution of the axis ratios with time is linear for almost all 
halos. Figure 18 shows the fractional rate of change, (b/a) /{b/a) and (c/a) /{c/a), of the 
axis ratios as a function of the value of the axis ratio in the final snapshot. The median and 
standard deviation of the distribution of {b/a) /{b/a) are 0.0089 h Gyr^^ and 0.0349 h Gyr^^ 
respectively. For c/a, they arc 0.0093 h Gyr~^ and 0.0297 h Gyr^^ respectively. Therefore, 
there is a weak tendency for undisturbed halos to become more spherical with time. Most 
halos require several Gyr before their flattening changes significantly; there are, however, 
a few outliers with quite significant changes in their axis ratios. Figure 18 demonstrates 
that there is no trend of {b/a) /{b/a) or {c/a) /{c/a) with the value of the axis ratio except 
for the outliers with very high (low) values of b/a, which could not have such high (low) 
rates of change if the values of b/a were not very high (low) in the final snapshot. We 
find no trend with any other halo property such as mass, spin parameter, pattern speed, 
substructure fraction, or alignment of the figure rotation axis with the angular momentum 
vector or minor axis. 



5. Conclusions 

We have detected rotation of the orientation of the major axis in most undisturbed halos 
of a ACDM cosmological simulation. The axis around which the figure rotates is very well 
aligned with the minor axis in most cases. It is also usually well aligned with the angular 
momentum vector. The distribution of pattern speeds is well fit by a log normal distribution, 

1 / 



with = 0-148 h km s"^ kpc"^ and a = 0.83. 

The pattern speed Vtp is correlated with spin parameter A. The median pattern speed 
rises from 0.12 h km s~^ kpc~^ for halos with A < 0.02 to 0.44 h km s~^ kpc~^ for halos 
with A > 0.06, with a spread of almost an order of magnitude around this median at a 
given value of A. The 11% of halos in the sample with the highest pattern speeds, Q,p > 
0.4 h km s~^ kpc~^, not only have large spin parameters, but also show particularly strong 
alignment between their figure rotation axes and their angular momentum vectors. There is 
no obvious correlation of the figure rotation properties with mass. The pattern speed and 
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figure rotation axis is coherent in the outer regions of the halo. Within 0.12 ryir, the pattern 
speed drops, particularly for those halos with fast figure rotation, and the internal alignment 
of the figure rotation axis deteriorates. This is probably an artifact of the numerical force 
softening. 

BFPM99 hypothesized that the spiral structure in NGC 2915 is due to figure rotation 
of a triaxial halo. The required pattern speed of 7 ± 1 km s~^ kpc^^ (MB03) is much higher 
than the pattern speeds seen in the simulated halos, and is estimated to have a probability 
of 5 X 10^^. We therefore conclude that the figure rotation of undisturbed ACDM halos is 
not able to produce this spiral structure. Halos with large values of A tend to have more 
substructure (Barnes & Efstathiou 1987), so there is a deficiency of halos with very high A in 
our sample. Because correlates with A, we cannot exclude the possibility that there exist 
halos with very high A whose figures rotate sufficiently quickly. However, halos with such 
high A are themselves very rare (MB 03), and if such halos fall out of our sample due to the 
presence of strong substructure, the effects of the substructure on the gas disk of NGC 2915 
would be of more concern than the slow rotation of the halo figure, a possibility BFPM99 
rule out due to the lack of any plausible companion in the vicinity. 

More generally, Bekki & Freeman (2002) found very weak if any enhancement of spiral 
structure in disk simulations with triaxial figures rotating at 0.77 km s~^ kpc~^, a value 
similar to the highest pattern speed seen in our sample. Therefore, it is unlikely that triaxial 
figure rotation can be detected by looking for spiral structure in extended gas disks. 

We have found that the axis ratios of undisturbed halos tend to become more spherical 
with time, with median fractional increases in the h/a and c/a axis ratios of 0.009 h Gyr"-^. 
The distributions of (6/ a) / (6/ a) and (c/ a) j [cj a) are relatively wide, with standard deviations 
of ~ 0.03 h Gyr~^. A few outliers have axis ratios that change quite significantly over the 
span of 1 Gyr. The rate of change of the axis ratios is not correlated with any other halo 
property. 
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Fig. 18. — Fractional rate of change of the h/a (left) and c/a (right) axis ratio of the halos 
in the sample, as a function of the z = Q axis ratio. 



